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Abstract 


The  Acoustic  Battlefield  Aid  (ABFA)  is  a  prototype  decision  aid  for 
assessing  the  performance  of  acoustic  sensors  in  different  environments. 
ABFA  combines  accurate  modeling  of  atmospheric  effects  on  acoustic 
signals  with  newly  developed  methods  for  determining  acoustic  array 
performance.  The  report  describes  the  methods  by  which  ABFA  calculates 
such  quantities  as  transmission  loss,  signal-to-noise  ratio,  detection 
probability,  array  direction-finding  accuracy,  and  position-finding 
accuracy  from  array  networks.  Some  example  calculations  are  provided 
for  propagation  over  rough  terrain  in  various  weather  conditions. 
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1.  Introduction 


Acoustical  systems  can  be  used  to  detect,  locate,  and  classify  targets  on  the 
battlefield.  They  offer  many  potential  benefits  for  such  uses:  low  cost, 
ruggedness,  versatility,  passive  operation,  and  detection  capability  in 
non-line-of-sight  situations.  They  also  provide  information  that  can 
complement  existing  optical  and  infrared  technologies.  Thus  it  appears 
that  acoustical  systems  will  play  a  valuable  role  on  the  future  battlefield. 

One  important  characteristic  of  acoustical  systems,  which  must  be 
accounted  for  in  battlefield  applications,  is  that  their  performance 
depends  heavily  on  atmospheric  conditions  and  the  local  terrain.  The 
atmosphere  and  terrain  alter  the  acoustic  signal  as  it  propagates  from  the 
source  to  the  sensor.  Some  important  propagation  effects  are  illustrated  in 
figure  1.  Unfortunately,  few  practical  tools  have  as  yet  been  available  that 
can  take  these  propagation  effects  into  account.  This  report  describes  a 
prototype  decision  aid,  called  the  Acoustic  Battlefield  Aid  (ABFA),  which 
has  been  designed  to  remedy  this  shortcoming. 

In  its  current  form,  ABFA  is  best  described  as  a  demonstration  model. 
Developed  by  the  U.S.  Army  Research  Laboratory  (ARL),  ABFA  provides 


Figure  1.  Influences  of 
atmosphere  and  terrain  on 
sound  propagation. 


Acoustic 


1 


a  new  conceptual  framework  for  how  an  acoustical  tactical  decision  aid 
can  be  designed  and  implemented.  Work  must  still  be  done  to  verify  its 
predictions,  and  to  refine  it  so  that  it  is  useful  in  Army  applications.  It  is 
hoped  that  the  general  design  and  underlying  methodologies  of  ABFA 
will  help  lead  future  development  of  useful,  flexible,  and  reliable  tools  for 
characterizing  acoustic  sensor  performance. 

ABFA  incorporates  recent  scientific  and  technological  advances  for 
predicting  acoustic  propagation,  and  for  characterizing  sensor  array 
performance  in  the  randomly  varying  (turbulent)  atmosphere.  Several 
state-of-the-art  research  products  are  merged  into  a  simple,  graphical  user 
interface.  The  acoustical  algorithms  have  been  carefully  implemented  to 
characterize  sensor  performance  in  a  computationally  efficient,  yet 
reliable,  maimer. 

The  emphasis  of  this  report  is  on  the  various  algorithms  used  for  ABFA's 
calculations.  It  is  not  intended  as  a  user's  guide,  nor  as  a  discussion  of  the 
performance  of  actual  sensor  systems.  These  matters  (although  of  course 
important)  are  beyond  the  scope  of  this  report. 

The  report  begins  with  a  discussion  of  some  of  the  broad  conceptual  issues 
relevant  to  the  design  of  ABFA  (sect.  2).  I  then  discuss  ABFA's  method  for 
handling  different  types  of  sources  and  receivers  (sect.  3),  atmospheric 
modeling  (sect.  4),  and  terrain  modeling  (sect.  5).  The  acoustic 
propagation  models  incorporated  into  ABFA  are  discussed  in  section  6, 
and  the  methods  by  which  sensor  performance  is  quantified  in  section  7. 
Section  8  provides  some  example  calculations. 

The  report  has  been  organized  so  that  readers  who  are  interested  mainly 
in  the  capabilities  of  ABFA,  rather  than  in  the  technical  details  of  its 
underlying  algorithms,  can  skip  directly  to  the  examples  in  section  8. 
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2.  The  General  Design  of  ABFA 


ABFA  is  ambitious  in  scope:  it  attempts  to  seamlessly  incorporate 
atmospheric  and  terrain  characterization,  the  properties  of  different 
acoustic  sources  and  receivers,  several  acoustic  propagation  models,  and 
statistical  methods  for  determining  the  effects  of  random  atmospheric 
fluctuations  (turbulence)  on  the  information  yielded  by  the  sensors.  In 
order  to  achieve  this  goal,  ABFA  makes  certain  reasonable  sacrifices 
regarding  the  exactness  of  the  calculations.  This  is  necessary  because  exact 
characterization  of  the  performance  of  a  given  sensor  requires  mimicking 
the  sensor's  electronics,  its  data  processing  algorithms,  and  the 
complicated  interactions  between  the  sensor  and  and  its  environment. 
Accounting  for  all  these  aspects  of  the  problem  would  be  a  prohibitively 
complex  and  computationally  infeasible  endeavor  for  a  multipurpose  tool 
such  as  ABFA.  Furthermore,  such  exactitude  would  involve  classified 
details  of  a  sensor's  design  and  operation.  Hence  there  is  a  strong 
motivation  to  decide  which  features  of  the  combined  acoustic 
propagation/ sensor  scenario  most  significantly  affect  the  sensor 
performance,  and  to  model  these  according  to  general  principles.  The 
following  are  some  of  the  general  modeling  principles  implicit  in  ABFA, 
and  their  corresponding  rationales: 

•  Principle:  Characterization  of  the  acoustic  sensor  performance  must 
account  for  the  influence  of  atmospheric  conditions  and  terrain  on 
signal  propagation. 

Rationale:  Sound  levels  are  dramatically  affected  by  refraction  and 
diffraction  by  atmospheric  wind  and  temperature  gradients,  and  by 
shadows  created  by  terrain.  Atmospheric  winds  also  produce  noise 
on  microphones  that  is  often  the  limiting  factor  in  low-frequency 
detection  capability.  For  propagation  distances  more  than  a  few 
hundred  meters,  acoustic  sensor  performance  depends  so  heavily  on 
atmospheric  conditions  and  terrain  that  a  decision  aid  not 
incorporating  them  is  nearly  useless. 

•  Principle:  Sensor  performance  should  be  characterized  by  metrics 
that  are  broadly  applicable  to  well-designed  sensors,  rather  than  by 
an  attempt  to  mimic  the  specific  electronics  and  data  processing 
algorithms  of  the  sensor. 
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Rationale:  Signal  processing  algorithms  used  by  acoustic  sensor 
systems  vary  widely  and  are  often  quite  complicated.  Furthermore, 
mimicking  most  good  processing  algorithms  would  impose  an 
impractical  computational  burden  on  the  propagation  model. 
Fortunately,  there  are  well-established  methods  for  approximating 
the  performance  obtainable  by  a  well-designed  sensor  and 
processing  algorithm.  The  metric  used  by  ABFA  is  called  the 
Cramer-Rao  lower  bound  (CRLB),  which  describes  the  performance  of 
an  optimal  sensor.  The  CRLB  is  fairly  easy  to  calculate,  and  previous 
studies  have  shown  that  good  processing  algorithms  attain 
performance  comparable  to  the  CRLB  (Song  and  Ritcey,  1996). 

•  Principle:  The  decision  aid  should  interface  simply  with  a  variety  of 
data  sources  for  the  atmospheric  and  terrain  inputs. 

Rationale:  The  decision  aid  may  be  used  in  scenarios  ranging  from 
simulations  to  actual  implementation  on  the  battlefield.  It  should 
interface  sensibly  with  available  data  sources,  regardless  of  the 
quality  and  resolution  of  the  input.  The  current  version  of  ABFA  is 
configured  to  accept  new  data  formats  without  changes  to  the 
functionality. 

•  Principle:  The  decision  aid  should  be  a  flexible,  multipurpose  tool, 
incorporating  calculations  and  displays  of  such  diverse  quantities  as 
transmission  loss,  signal-to-noise  ratio  (SNR),  detection  probabilities, 
accuracy  of  sensor  arrays  for  determining  target  bearings,  and 
accuracy  of  multiple  array  configurations  for  locating  targets. 

Rationale:  The  most  difficult  part  of  the  acoustical  calculation  is  the 
effect  of  atmospheric  conditions  and  terrain  on  the  acoustic 
propagation.  Other  types  of  plots  are  basically  enhancements  of  the 
propagation  information,  with  some  small  additional  processing. 
Therefore  there  is  little  reason  for  making  separate  decision  aids  for 
different  display  purposes. 

With  regard  to  the  last  principle,  the  main  effect  of  the  weather  and  terrain 
is  to  alter  the  transmission  loss  (also  called  the  propagation  loss)  from  the 
source  to  the  receiver.  Calculation  of  the  transmission  loss  is  the  first  step 
in  all  ABFA  calculations.  Once  the  transmission  loss  is  known,  the  SNR 
can  be  calculated  from  the  source  and  background  noise  levels.  Other 
calculations,  involving  array  performance,  follow  from  the  SNR.  The  basic 
hierarchy  in  calculations  made  by  ABFA,  and  the  inputs  at  each  step,  are 
shown  in  figure  2.  Results  of  any  of  the  intermediate  calculations  are 
available  for  plotting. 
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ABFA  was  programmed  in  the  high-level  Matlab  language.  The  built-in 
matrix  routines,  graphical  capabilities,  and  GUI  (graphical  user  interface) 
functions  of  Matlab  made  it  possible  to  produce  ABFA  with  limited 
programming  effort.  The  initial  screen  presented  to  the  user  is  shown  in 
figure  3. 


Figure  2.  Hierarchy  of 
calculations  available 
from  ABFA.  Inputs: 

G  =  ground  and  terrain; 
N  =  background  noise; 
S  =  sensor  (receiver); 

T  =  target  (source);  and 
W  =  weather 
(atmospheric 
characterization) . 
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Figure  3.  ABFA's  graphical 
user  interface.  This  initial 
screen  shows  topography 
and  all  source  ( x )  and 
receiver  (Q)  locations. 
Alternative  displays,  such 
as  signal-to-noise  ratio  and 
probability  of  detection, 
are  selected  from 
"Display"  menu. 
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The  drawback  of  using  Matlab  to  create  ABFA  is,  of  course,  that  the 
Matlab  interpreter  must  be  available  to  run  the  program.  Current  ART 
projects  include  the  creation  of  an  attractive  GUI,  written  in  the  C++ 
language,  for  rurming  acoustical  prediction  models.  Plans  include 
conversion  of  many  of  ABFA's  algorithms  to  C++.  With  regard  to  these 
plans,  it  is  worth  pointing  out  that  The  MathWorks,  Inc.  (makers  of 
Matlab),  recently  introduced  a  compiler  that  converts  Matlab  code  to  C. 
However,  routines  with  graphics  commands  cannot  be  compiled. 
Therefore  only  the  computational  parts  of  the  code  can  be  converted  to  C 
with  the  Matlab  compiler. 
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3.  Source  and  Receiver  Models 


As  discussed  above,  one  of  the  guiding  principles  in  the  design  of  ABFA 
was  that  it  should  estimate  the  performance  of  a  well-designed  sensor,  rather 
than  predicting  precisely  how  a  particular  sensor  system  will  operate.  The 
decision  to  design  ABFA  in  this  manner  is,  of  course,  a  compromise:  it 
would  be  ideal  if  the  performance  of  a  given  sensor  could  be  simulated 
precisely,  but  practical  considerations  often  prevent  this. 

In  ABFA,  each  source  is  assumed  to  be  omnidirectional  and  is  represented 
by  a  single  characteristic  frequency  fc  and  a  finite  bandwidth  A/.  The 
characteristic  frequency  is  taken  to  be  the  geometric  mean  of  the  lower  {fi) 
and  upper  {fu)  frequencies  in  the  band: 

fc  =  Vfefu, 

=  fu  -  fe- 

To  calculate  the  signal  transmission  loss  (the  difference  between  the  signal 
amplitude  at  a  distance  of  1  m  from  the  source,*  and  its  level  at  the 
receiver),  the  frequency  fc  is  used.  The  source  is  essentially  modeled  as  a 
harmonic  point  source  for  this  part  of  the  computation.  But  for  calculating 
the  SNR,  all  the  noise  energy  in  the  band  [fe,  fu]  is  summed  to  produce  the 
total  noise  background. 

The  reason  for  using  a  single  frequency  and  bandwidth  to  model  sources 
is  to  save  computational  time.  Most  of  the  computational  time  in  an  ABFA 
calculation  is  usually  spent  in  the  propagation  model,  which  must  be  run 
separately  for  each  frequency.  Hence  the  single-frequency  method 
minimizes  computational  time. 

It  should  be  made  clear  that,  even  though  ABFA  runs  the  propagation 
model  at  just  one  frequency,  there  is  no  implication  that  the  actual  source 
emits  acoustic  waves  just  at  that  single  frequency.  What  is  implied  is  that 
the  propagation  characteristics  of  that  single  frequency  are  representative 
of  propagation  for  the  frequency  band  (or  set  of  frequencies)  used  by  a 
particular  sensor.  For  example,  suppose  a  tank  is  modeled  as  a  source 
having  center  frequency  fc  =  100  Hz  and  bandwidth  A/  =  150  Hz  {fe  = 

50  Hz,  fu  =  200  Hz).  Further  suppose  that  a  particular  sensor  and 

*The  reference  distance  of  1  m  from  the  source  is  an  accepted  standard  in  acoustics 
research  and  noise  control  applications. 
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processing  algorithm  uses  information  from  frequency  lines  such  as  80, 
100,  and  150  Hz.  Since  the  transmission  losses  at  these  three  frequencies 
are  normally  similar  to  those  at  100  Hz,  the  ABFA  prediction  should  work 
reasonably  well.  On  the  other  hand,  if  information  at  a  very  low  frequency 
(such  as  20  Hz)  is  used  by  the  sensor,  the  ABFA  propagation  calculation 
may  not  be  satisfactorily  representative.  Furthermore,  contributions  from 
wind  noise  below  =  50  Hz  would  be  neglected,  which  may  be  quite 
significant  in  high-wind  conditions. 

The  single-frequency  method  used  by  ABFA  is  just  a  convenient  (and 
usually  reasonable)  assumption  used  to  create  the  prototype.  Eventually  it 
may  be  desirable  to  upgrade  ABFA  from  its  present  form  to  one  that 
models  the  performance  of  specific  sensor  systems  with  higher  accuracy, 
in  which  case  the  single  representative  frequency  method  can  easily  be 
upgraded  to  incorporate  several  representative  frequencies. 

In  the  current  version  of  ABFA,  receivers  consist  of  ideal  point 
microphones  having  flat  (frequency-independent)  responses.  For  most 
microphone  types  and  acoustic  frequencies  of  interest,  the  assumptions  of 
a  point  receiver  and  a  flat  response  are  quite  reasonable.  The  microphones 
may  be  arranged  in  coherent-processing  arrays  having  a  linear,  circular,  or 
triangular  geometry.  In  all  these  cases,  the  number  of  elements  and  the 
spacing  are  chosen  by  the  user. 

The  microphones  are  assumed  to  be  equipped  with  9-cm-diameter 
windscreens.  The  effect  of  the  windscreen  is  accounted  for  in  the  SNR 
calculation  (sect.  7.1).  Arbitrarily  sized  windscreens  may  be  incorporated 
into  a  future  version  of  ABFA. 
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4.  Atmospheric  Models 


When  the  wind  and  temperature  vary  with  height,  sound  energy  is 
refracted,  a  process  that  tremendously  influences  sound  levels  in  the 
atmosphere.  ABFA  provides  three  basic  methods  for  entering  the  mean 
wind  and  temperature  profiles  used  for  the  calculations.  The  first  of  these 
involves  user  selection  from  a  half-dozen  "typical"  cases,  described  in 
detail  by  Noble  (1994).  The  second  method,  called  Monin-Obukhov  (MO) 
similarity  scaling,  is  based  on  a  statistical  parameterization  for  atmospheric 
turbulence  that  has  been  well  verified  during  the  past  several  decades. 
MO-similarity  is  described  in  section  4.1.  The  third  method  is  simply  to 
load  a  file  containing  arbitrary  profiles.  This  method  is  useful  if  data  are 
available  from  a  radiosonde  or  other  meteorological  instrumentation. 

ABFA's  calculations  involving  arrays  require  statistical  information  on  the 
turbulent  fluctuations  of  wind  and  temperature,  in  addition  to  the  plain 
mean  profiles.  The  turbulent  fluctuations  are  parameterized  with  structure 
functions,  as  described  in  section  4.2. 


4.1  Monin-Obukhov  Similarity  Modeling 

The  wind  and  temperature  profiles  in  the  atmospheric  surface  layer  (the 
10  to  200  m  nearest  the  ground)  are  determined  by  turbulent  mixing. 
Although  it  is  not  generally  possible  to  solve  problems  involving 
turbulent  mixing  on  the  basis  of  first  principles,  there  are  scaling  methods 
that  have  been  found  to  yield  reliable  results.  For  the  atmospheric  surface 
layer  profiles,  the  MO-similarity  (Panofsky  and  Dutton,  1984;  Stull,  1988) 
scaling  method  is  often  used. 

MO-similarity  is  based  on  the  premise  that  atmospheric  turbulence 
structure  near  the  ground  depends  entirely  on  four  quantities:  the  height 
z,  the  stress  r  exerted  by  the  wind  on  the  ground,  the  surface  heat  flux  Q, 
and  the  buoyancy  (3  =  glOo,  where  g  is  gravitational  acceleration  and  6q  is 
the  surface  temperature.  From  these  four  quantities,  representative  scales 
for  length,  velocity,  and  temperature  can  be  formulated. 

The  velocity  scale,  called  the  friction  velocity,  is 


rt*  = 


(1) 
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(2) 


where  p  is  the  density  of  air.  The  surface-layer  temperature  scale  is 

rr  -  Q 
1^  — 

Last,  the  length  scale,  called  the  Monin-Obukhov  length,  is 


Lm,o  — 


Wi 


Kj3Q 


(3) 


where  n  =  0.4  is  called  von  Karman's  constant. 


According  to  MO-similarity,  the  mean  profiles  must  be  universal  functions 
of  the  ratio  C  =  zjLmo-  As  a  result,  the  profiles  for  the  mean  wind  speed 
and  temperature  can  be  written 


u{z) 

T{z) 


K 


In  ^ —  T'm 
L  zo 


+  'I'm 


Lr 


and 


T  (Zref)  + 


T„Pt 

K 


^ref 

Lmo 


(4) 

(5) 


where  zq  is  called  the  roughness  height  (height  where  the  mean  wind 
vanishes),  Zref  is  some  reference  height  at  which  the  temperature  is 
known,  Pt  ~  0.74,  'kjvf  is  the  universal  profile  function  for  momentum, 
and  'I'h  is  the  universal  profile  function  for  heat.  For  the  universal 
functions,  ABFA  uses  the  forms  proposed  by  Carl  et  al  (1973;  see  also 
L'Esperance  et  al,  1993): 

4'm(C)  =  4'h(C)  =  ^  - h  \/3tan“^  C  <  0, 

'Em  (C)  =  'I'H  (C) = -4.7C,  C>0, 


where 

fMiC)  =  (1-15C)-'/", 

=  (l-9C)-'/^ 


Sometimes,  particularly  when  conditions  are  very  sunny  and  the  wind 
very  light,  MO-similarity  is  known  to  break  down.  The  reason  is  that 
large,  boundary-layer-size  eddies  (typically  1  to  2  km  in  size)  begin  to  play 
an  important  role  in  turbulent  mixing  even  very  near  the  ground.  In  these 
cases,  instead  of  using  the  scaling  parameters  described  above,  one  should 
switch  to  the  convective  boundary-layer  scaling  parameters  (Stull,  1988) 


0^ 
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These  scales  are  for  velocity  and  temperature,  respectively  The 
appropriate  length  scale  is  Zi,  the  height  of  the  boundary-layer  inversion. 
In  situations  where  turbulent  structure  depends  primarily  on  free 
convection,  ABFA  automatically  switches  from  MO-similarity  to 
convective  scaling. 

ABFA  allows  the  user  to  enter  arbitrary  values  for  u^,  T*,  and  z*. 
Alternatively,  the  user  can  select  a  set  of  values  based  on  a  subjective 
description  of  the  weather  conditions.  Table  1  shows  the  different  weather 
conditions  available,  and  the  corresponding  values  of  u*,  T*,  and  Zi. 

4.2  Structure-Function  Parameters 

Structure-function  parameters  are  useful  for  characterizing  the  "strength" 
of  turbulence.  They  appear  in  statistical  descriptions  of  wave  propagation 
when  scattering  by  eddies  in  the  inertial  subrange  is  analyzed.  (The 
inertial  subrange  in  the  atmosphere  generally  consists  of  eddies  smaller 
than  the  height  from  the  ground,  and  larger  than  1  mm.)  The 
structure-function  parameters  are  given  by  the  equations 

Cf.  = 

C\  = 

(as  given  in  Stull,  1998,  and  elsewhere),  where  the  dissipation  rate  for 
turbulent  kinetic  energy  is 


o  Q  N 

ut  w'r 


€  ~  max 


y0.4z’  Zi  j  ’ 

and  the  destruction  rate  for  temperature  variance  is 


X  ~  0.43- 


Table  1.  Preselected 
meteorological  conditions,  - 
based  on  Monin-Obukhov 

Description 

u*  (m/s) 

A(K) 

2i  (m) 

mostly  sunny  and  calm 

0.1 

-0.3 

1000 

similarity  scaling,  available 
in  ABFA. 

mostly  sunny  and  windy 

0.5 

-0.1 

1000 

clear  night  and  calm 

0.1 

0.02 

100 

clear  night  and  windy 

0.5 

0.05 

100 

overcast  (day  or  night)  and  calm 

0.1 

-0.01 

1000 

overcast  (day  or  night)  and  windy 

0.5 

-0.01 

1000 
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5.  Terrain  Models 


ABFA  provides  two  basic  options  for  generating  topography  (elevation 
models)  in  the  simulation.  One  of  these  is  simply  to  load  a  file  that 
contains  a  two-dimensional  array  whose  elements  correspond  to  heights 
at  evenly  spaced,  rectangular  grid  points.  For  the  other  option — 
generating  a  random  terrain  based  on  a  height  spectrum — ABFA  provides 
two  spectra:  Gaussian  and  von  Karman.  The  method  of  generating  the 
synthetic  terrain  is  essentially  the  same  for  either  model.  First,  a  power 
spectrum  for  the  heights  is  used  to  determine  the  strength  of  the  Fourier 
(wavenumber-domain)  modes.  Each  Fourier  component  is  then 
multiplied  by  a  random  phase.  Finally,  an  inverse  Fourier  transform  is 
applied  to  produce  the  spatial-domain  terrain.  This  section  describes  the 
two  prescribed  spectra,  along  with  the  selections  available  for  small-scale 
terrain  roughness  features,  such  as  grass  and  trees. 


5.1  Gaussian  Spectrum 

The  Gaussian  model  uses  the  following  power  spectrum: 


S{k) 


—  exp 


(6) 


where  k  is  the  radial  wavenumber  (/c^  =  k‘^  +  ky),  is  the  variance  in  the 
terrain  height,  and  L  is  the  length  scale  of  the  terrain  features.  The 
Gaussian  model  yields  terrain  having  a  smooth  appearance,  with  the 
"hills"  all  having  horizontal  dimensions  comparable  to  L. 


5.2  Von  Karman  Spectrum 

The  von  Karman  model  is  based  on  the  following  equation  for  the  power 
spectrum,  which  is  similar  to  one  commonly  used  to  describe  turbulence: 

=  (7) 


12 


Here  the  parameter  v  is  called  the  spectral  roll-off  parameter.  Like  the 
Gaussian  model,  the  von  Karman  model  produces  terrain  with  strong 
features  having  length  scales  comparable  to  L.  The  main  difference 


between  the  models  is  that  the  von  Karman  model  also  contains  noticeable 
features  having  scales  smaller  than  L,  resulting  in  a  realistic,  rough 
appearance  to  the  terrain  not  evident  in  the  Gaussian  model.  The  amount 
of  roughness  can  be  controlled  by  adjustments  to  v,  which  is  related  to  the 
fractal  dimension  of  the  terrain.  Higher  values  of  v  correspond  to  more 
roughness  and  higher  fractal  dimension. 

5.3  Surface  Roughness  Characterization 

Besides  the  actual  changes  in  terrain  elevation,  it  is  important  to  have 
some  characterization  of  the  ground  surface  and  the  small-scale 
aerodynamic  features  existing  on  the  terrain  (e.g.,  desert,  forest,  urban).  In 
the  terminology  of  fluid  dynamics,  these  features  are  called  the  roughness 
elements.  They  determine  the  roughness  length  zq  used  in  the  wind  profile, 
equation  (4).  ABFA  provides  a  set  of  choices  for  the  terrain  ground  type, 
each  having  a  characteristic  roughness.  The  terrain  types  and 
corresponding  roughness  lengths  are  shown  in  table  2.  The  roughness 
lengths  were  adapted  from  Panofsky  and  Dutton's  (1984)  table  6.2. 


Table  2.  Selections  for 
terrain  type  and 
corresponding  values  for 
surface  roughness  length 

Description 

Roughness  length 

20  (m) 

asphalt  or  ice 

10-^ 

zo- 

open  water 

10“® 

snow 

2  X  10“® 

desert  (sand) 

10“^ 

short  grass 

10“^ 

long  grass  or  crops 

5  X  10“^ 

brush 

0.1 

suburban 

0.4 

forest 

1 

city 

10 
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6.  Propagation  Models 


Calculations  of  sensor  array  performance  depend  on  knowledge  of  the 
signal-to-noise  ratio  (SNR)  and  the  signal  coherence  between  the  sensors 
in  the  array.  In  order  to  find  the  SNR,  one  needs  to  know  the  reduction  in 
signal  level  due  to  propagation  between  the  source  and  the  receiver.  When 
expressed  in  decibels,  this  reduction  is  called  the  transmission  loss  (TL).  It  is 
defined  by 

TL  =  201og5-201ogio|p(R)|,  (8) 

where  p  (R)  is  the  complex  acoustic  pressure  at  the  location  R,  and  S  is  the 
strength  of  the  source  measured  in  free  space  at  1  m.  (By  definition,  the  TL 
of  a  source  radiating  into  free  space,  measured  at  a  distance  1  m,  is  0  dB.) 

ABFA  contains  five  models  (four  distinct  and  one  hybrid)  for  determining 
TL.  These  are  discussed  in  sections  6.1  to  6.5,  in  order  of  increasing 
complexity.  The  signal  coherence  between  the  sensors  is  quantified  via  the 
mutual  coherence  function  (MCF),  which  is  discussed  in  section  6.6. 

6.1  Hard  Surface  Model 

The  sound  field  above  a  perfectly  reflecting  (hard)  surface  can  be  written 
as  the  sum  of  a  contribution  traveling  directly  from  the  source  and  one 
reflecting  from  the  surface.  The  surface  reflection  can  be  thought  of  as 
emanating  from  an  image  source  at  a  height  —Zg,  where  Zg  is  the  height  of 
the  actual  source,  as  shown  in  figure  4.  The  acoustic  pressure  field 
generated  by  a  single,  unit-strength  source  is 

(9) 

K 

where  k  is  the  wavenumber,  and  R  is  the  distance  between  the  source  and 
receiver.  The  function  G  (R)  is  called  the  free-space  Green's  function.  When 
the  contributions  from  the  real  and  image  sources  are  summed  and  then 
multiplied  by  the  source  strength,  we  arrive  at  the  following  equation  for 
the  total  pressure  field: 

p  =  S[G{R)  +  G{R,)],  (10) 

where 

R?  =  {x-  Xgf  +  {y-  ygf'  +  {z  -  Zg)"^ , 

Rf  =  {x-  Xgf  +  (y  -  ygf  +  {z  +  Zgf  . 
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Figure  4.  A  source  above  a 

hard  boundary,  showing  actual  source 

resulting  image  source. 


receiver 


{x,  y,  z) 


R 


air 


ground 


image  source 

{xg,  ys,  -Zs) 


6.2  Impedance  Ground  Model 


The  hard-surface  model  discussed  in  the  previous  section  assumes  that  all 
the  sound  energy  incident  on  the  ground  plane  is  reflected.  In  actuality, 
the  ground  absorbs  some  fraction  of  the  sound  energy.  We  can  account  for 
this  partial  absorption  by  modeling  the  ground  as  a  finite  impedance 
surface.  A  perfectly  hard  surface  corresponds  to  an  infinite  impedance. 

Chien  and  Soroka  (1975)  have  previously  determined  the  sound  field 
above  a  locally  reacting,  finite-impedance  surface.  Their  solution  has  a 
form  very  similar  to  the  hard-surface  case: 


p  =  S[G{R)  +  Q{i;,w)G{Ri)] 


(11) 


where  Q  {ip,  w)  modifies  the  strength  of  the  reflection.  The  angle  ip  is 
defined  in  figure  4,  and  the  quantity  w  is  given  by 


(12) 


where  po  is  the  air  density,  cq  the  sound  speed  in  the  air,  and  Zg  the  surface 
impedance.  The  function  Q  {ip,  w)  is  written 


Q  {ip,  w)  =  R{ip)  +  [1  -  R  {ip)]  F  {w) 


(13) 


in  which 


sin-0  -  pqCqIZs 
sin-f/;  -F  poco/Zg 


(14) 


is  called  the  plane-wave  reflection  coefficient,  and 


(15) 
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ABFA  determines  the  surface  impedance  using  a  Biot/ Allard-like,  porous, 
relaxational  model  (Wilson,  1997b).  There  are  four  parameters  in  the 
model:  the  static  flow  resistivity  a,  the  tortuosity  q,  the  porosity  fl,  and  the 
Biot/ Allard  shape  factor  sb-  Values  for  these  parameters  are  set  based  on 
the  terrain  ground  type  (sect.  5.3),  and  are  given  in  table  3.  These  values 
are  representative  of  the  measurements  described  by  Embleton  et  al. 

(1983),  Martens  et  al.  (1985),  and  Attenborough  and  Buser  (1988). 


Table  3.  Ground 
parameters  used  in 
relaxational  model.  Water 
and  ice  cases  both  use  an 
infinite  surface  impedance. 


Description 

(T  (Pa  •  s/m^) 

<1^ 

n  (%) 

SB 

asphalt,  city 

3  X  10'^ 

3.2 

10 

1 

open  water  or  ice 

infinite 

— 

— 

— 

snow 

10‘‘ 

1.7 

60 

1 

desert  (sand) 

4  X  10® 

1.6 

40 

1 

grass  (short  or  long),  crops,  brush,  suburban 

2  X  10® 

1.4 

50 

1 

forest 

10® 

1.3 

60 

1 

6.3  Wedge  Diffraction  Model 

The  propagation  models  discussed  so  far  have  neglected  terrain  features 
such  as  hills  and  valleys.  Although  there  are  no  exact,  analytical  solutions 
for  the  sound  field  over  arbitrary  terrain  features,  solutions  have  been 
calculated  for  some  idealized  surface  geometries.  One  such  geometry  that 
has  a  tractable  solution  is  the  wedge  shape,  shown  in  figure  5.  The  angle  (3 
in  figure  5  is  called  the  exterior  wedge  angle,  and  v  =  n/ (5  {>1/2)  is  called 
the  wedge  index. 


Figure  5.  Geometry  for 
propagation  around  a 
wedge.  Shown  is  a  cross 
section  through  wedge;  it 
is  actually  infinite  along 
axis  perpendicular  to  page 
(z-axis). 
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PD  = 


The  form  of  the  solution  for  the  wave  field  around  a  wedge  that  is 
incorporated  into  ABFA  is  due  to  Pierce  (1981).  The  total  pressure  field  is 
decomposed  into  two  contributions,  one  part  from  geometrical  acoustics 
and  the  other  from  diffraction: 


P  =  PGA  +  PD- 


One  can  obtain  Pierce's  form  of  the  geometrical  acoustics  solution  by 
writing  the  Green's  function  for  a  point  source  in  a  slightly  different  form. 
First  note  from  figure  5  and  the  law  of  cosines  that 

+  Tg  —  2rrs  cos  (  +  {z  —  Zg)'^  , 


where  (  =  (p  —  fg  the  angle  between  the  source  and  receiver.  Let  us  now 
write  the  Green's  function  in  terms  of  dependence  on  (  instead  of 
dependence  on  R: 

exp  i ik\J -\-  rl  —  2rrs  cos  C  +  (z  —  | 

g{C)  =  G{R)  =  — \  ^  (16) 

_|_  j,2  _  2rrs  cos  C  +  (2:  —  Zg)‘^ 

Using  this  form  of  the  Green's  function,  we  find  that  the  geometric 
acoustics  component  of  the  solution  is  (Pierce's  equation  (9-8.1)) 


PGA 


S 


'^Q  {2(5n- 4>  + (pg) +  ^g  {2j3n- (p  -  (pg)  , 


(17) 


where  the  summation  extends  over  all  values  for  which  the  argument  of 
the  Green's  function  is  between  — vr  and  tt.  It  is  straightforward  to  show 
that  equation  (17)  reduces  to  (10)  when  f3  =  n. 

The  equation  used  by  ABFA  for  the  diffracted  contribution  is  an 
asymptotic  expression  valid  for  A:r  »  1  and  kvg  »  1,  given  by  Pierce  as 
equations  (9-8.10)  and  (9-8.11): 


5 


exp  {ikL  —  ITT  I  A)  sin  utt 
^/2L 


-  V - VT- 

1  —  cos  VTT  COS  V 

T-cjyy  1  - -J-/  V* - \'r  ‘  'T  ^ ; 

{(p  —  (pg)  1  —  cos  VTT  cos  u  {(p  +  (ps)  \  ’ 

(18) 

Ad{X)  = 

sign(X)  [f{\X\)-ig{\X\)], 

(19) 

M,{(P)  = 

cos  VTT  —  cos  ucp 

(20) 

U  {1  —  cos  UTT  cos  UCpf^'^ 

T^  = 

krvg  , 

— — ,  and 
ttL 

(21) 

= 

(r  -h  r*)^  -h  2^. 

(22) 
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where 


The  functions  /  and  g  are  called  the  Fresnel  auxiliary  functions,  and  are 
given  approximately  by  (Abramowitz  and  Stegun,  1965,  eq  (7.3.32)  and 
(7.3.33)) 

1  +  0.926X 

- and 

2  +  1.792a; +  3.104x2’ 

1 

2  +  4.142X  +  3.492x2  +  6.670x3  ' 


fix)  = 
gix)  = 


ABFA  provides  a  propagation  model  option  called  the  "wedge  model." 
The  wedge  model,  as  the  name  suggests,  is  based  on  the  solution  for  the 
field  around  a  wedge  given  above.  Given  source  and  receiver  locations, 
the  intervening  topographical  features  are  approximated  by  a  single 
wedge.  The  point  of  the  wedge  corresponds  to  the  high  point  (apex)  of  the 
intervening  terrain,  as  shown  in  figure  6.  Note  that  the  propagation 
coordinates  used  by  ABFA  (x  and  y  horizontal,  and  2;  vertical)  are  different 
from  the  coordinate  system  used  in  figure  5.  It  is  a  simple  matter,  however, 
to  shift  the  propagation  coordinates  so  that  the  point  of  the  wedge 
corresponds  to  the  origin,  and  then  rotate  them  to  agree  with  figure  5. 

6.4  Fast  Field  Program 

The  fast  field  program  (FFP)  is  a  method  for  computing  sound  fields  in 
horizontally  stratified  environments,  and  hence  can  be  used  to  calculate 
refraction  by  vertical  gradients  of  wind  and  temperature.  Originally 
developed  for  ocean  acoustics,  the  FFP  was  first  used  for  atmospheric 
propagation  over  a  decade  ago  (Raspet  et  al,  1985).  The  FFP  formulation 
incorporated  into  ABFA  is  described  elsewhere  (Wilson,  1993). 

The  FFP  is  an  excellent  tool  for  determining  meteorological  effects  on 
sound  propagation.  Unlike  such  methods  as  ray  tracing,  the  FFP  works 
well  at  low  frequency  and  in  refractive  shadow  regions.  Unfortunately, 
since  it  is  limited  to  horizontally  stratified  environments,  the  FFP  is 
applicable  only  to  propagation  over  flat  ground.  Therefore  the  FFP  can 
give  good  results  only  in  situations  where  the  terrain  is  nearly  flat. 


Figure  6.  Approximation 
of  rough  terrain  by  a 
wedge  shape. 


apex 
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The  FTP  in  ABFA  uses  the  meteorological  models  discussed  earlier 
(sect.  4).  The  program  automatically  selects  the  resolution  of  the 
calculation  in  a  manner  that  is  appropriate  for  the  dimensions  of  the 
simulated  grid.  The  ground  is  modeled  as  a  porous  half  space,  with 
parameters  given  by  table  3  (p  16). 

6.5  Hybrid  Wedge  Models 

In  order  to  handle  situations  where  both  terrain  and  atmospheric  effects 
are  important,  ABFA  contains  a  hybrid  propagation  model  that  uses  both 
the  wedge  and  FFP  calculations.  The  hybrid  model  is  reasonably  fast,  but 
not  as  rigorous  as  the  pure  wedge  or  FFP  models.  The  hybrid  model 
calculates  the  transmission  losses  corresponding  to  the  hard,  flat  surface 
(TL//),  the  wedge  (TLw)/  and  the  FFP  (TLppp),  and  then  combines  them 
incoherently: 

TLtotoi  =  TLiTirp  +  TLw  —  TLp.  (23) 

Note  that  the  transmission  losses  for  the  wedge  and  FFP  are  added 
together,  and  then  the  TL  for  the  hard  surface  is  subtracted.  The  rationale 
behind  this  approach  is  that  both  the  wedge  and  FFP  calculations  already 
include  the  effect  of  spherical  spreading  over  a  ground  surface.  Therefore 
the  hard-surface  contribution  from  the  wedge  model  needs  to  be 
subtracted  from  the  overall  solution.  The  remaining  propagation  effects 
are  then,  in  essence,  simply  summed  together  in  decibels.  These  effects 
include  terrain  (from  the  wedge  model),  meteorological  effects  (from  the 
FFP),  ground  impedance  effects  (from  the  FFP),  and  atmospheric 
absorption  (from  the  FFP). 

ABFA  also  incorporates  a  hybrid  wedge /impedance  ground  model  option 
for  convenience.  The  equation  for  this  model  is  the  same  as  equation  (23), 
although  with  TLg/  (the  TL  from  the  impedance  ground  model)  replacing 
TLppp.  The  hybrid  wedge /impedance  ground  model  is  reasonably  fast, 
but  neglects  atmospheric  effects. 

6.6  Mutual  Coherence  Function 

The  method  for  determining  coherence  between  signals  at  the  sensors  in  a 
receiving  array  is  rather  crude  in  comparison  to  the  methods  for 
determining  mean  levels  described  in  the  previous  sections.  The  main 
reason  for  this  shortcoming  is  a  lack  of  knowledge  of  atmospheric 
propagation  effects  on  second-order  moments  of  acoustic  fields.  Second 
moments  simply  have  not  received  as  much  research  attention  as  the  first 
moments  (means). 


19 


The  coherence  equation  used  by  ABFA,  which  is  valid  for  line-of-sight 
propagation  through  inertial-subrange  turbulence,  is  (Wilson,  1998b) 


r  (p)  =  exp 


7rk‘^R 

4 


0.464^  -h  3.40 


(24) 


Here  k  is  the  wavenumber,  R  the  total  propagation  distance,  Tq  the 
ambient  temperature  in  kelvins,  cq  the  ambient  sound  speed,  the 
average  structure-function  parameter  for  temperature,  Cy  the  average 
structure-function  parameter  for  velocity,  and  p  the  separation  between 
the  sensors.  The  average  structure-function  parameters  are  found  by 
integration  over  the  propagation  path: 

~  H  Jo 
^  H  Jo 

Height-dependent  expressions  for  (z)  and  Cy  (z)  are  given  in 
section  4.2. 
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7.  Sensor  Performance  Predictions 


7.1  Signal-to-Noise  Ratio 

The  SNR  provides  a  very  good  indication  of  sensor  performance. 
Generally,  if  the  SNR  is  well  above  0  dB,  the  sensor  will  function 
satisfactorily  On  the  other  hand,  negative  SNR  normally  indicates  that  the 
sensor  will  yield  little  useful  information.  The  SNR  (in  decibels)  simply 
equals  the  received  signal  level  (RL)  less  the  noise  level  (NL): 

SNR  =  RL  -  NL. 

The  received  signal  level  is  simply  the  sum  of  the  source  level  (SL)  (as 
would  be  measured  at  1  m  in  free  space)  and  the  transmission  loss  as 
determined  by  the  propagation  model: 

RL  =  SL  +  TL.  (25) 


The  most  important  noise  source  in  most  low-frequency  detection 
scenarios  is  wind  noise.  ABFA  models  the  spectrum  of  the  wind  noise 
(WL)  using  the  following  empirical  equation: 


WL  =  max  40  1  — 


log  /  -  log  20 
log  500  —  log  20 


+  2.7U,  0 


where  U  is  the  wind  speed  in  meters  per  second  (calculated  with  the 
atmospheric  model).  The  units  of  WL  are  given  in  decibels  re  10“^^  W  per 
one-third  octave  band.  Equation  (26)  was  determined  by  curve-fitting  to 
the  data  plotted  in  figure  5.60  of  Beranek  (1988).  These  data  were  recorded 
by  a  microphone  with  a  9-cm-diameter  windscreen.  The  following 
equation  relates  WL  in  decibels  to  the  corresponding  spectrum  Sw  if)  in 
power  units: 

WL  =  lOlogSwif)  ■ 


Besides  the  wind  noise,  ABFA  also  incorporates  an  environmental  noise 
background,  (/).  Various  values  are  assigned  to  Se  if)  depending  on 
the  user's  selection  from  such  choices  as  "rural  or  forest"  and  "light 


battle." 


To  find  the  total  noise  level,  ABFA  integrates  S'e  (/)  and  Sw  (/)  over  the 
bandwidth  of  the  source: 

NL  =  10  log  [Se  if)  +  Sw  (/)]  df.  (27) 

Jfi 
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7.2  Probability  of  Detection:  Single  Receiver 

In  most  situations,  the  SNR  at  a  receiver  undergoes  significant 
fluctuations.  These  fluctuations  can  be  attributed,  in  part,  to  random 
changes  in  the  background  noise  level,  or  to  changes  in  the  source  (target) 
power  output.  The  SNR  fluctuations  also  occur  as  a  result  of  random 
propagation  effects,  such  as  scattering  by  atmospheric  turbulence.  Because 
of  the  random  nature  of  the  SNR,  we  cannot  say  with  certainty  whether  a 
source  will  be  detected  in  a  given  environment.  Rather,  only  a  -probability  of 
detection  can  be  specified.  Actually,  there  are  four  possible  events  in  a 
detection  scenario,  each  having  its  own  associated  probability: 

1.  No  source  (target)  is  present,  but  the  detection  algorithm  incorrectly 
deduces  that  a  source  is  present.  This  is  called  a  false  alarm,  the 
probability  of  which  is  denoted  as  Pfa- 

2.  No  source  is  present,  and  the  detection  algorithm  correctly  deduces 
that  no  source  is  present.  The  probability  of  this  event  is  1  —  Pfa- 

3.  A  source  is  present,  and  the  detection  algorithm  correctly  deduces 
that  the  source  is  present.  This  is  called  a  detection,  with  probability 
Pd- 

4.  A  source  is  present,  but  the  detection  algorithm  incorrectly  deduces 
that  no  source  is  present.  This  is  called  a  false  dismissal,  with 
probability  Pfd  =  1  —  Pd- 

Obviously,  we  would  like  to  have  as  high  a  detection  probability  as 
possible,  while  undergoing  few  false  alarms.  These  two  desires  must  be 
balanced  against  one  another.  If  we  set  the  threshold  of  our  detector  to  a 
very  low  value,  there  will  be  a  lot  of  false  alarms,  but  we  have  an  excellent 
chance  of  detecting  a  real  target  when  one  comes  along.  A  high  threshold 
produces  fewer  false  alarms,  but  may  cause  us  to  miss  a  target.  A  popular 
method  for  formulating  this  trade-off,  and  the  method  used  by  ABFA,  is 
the  Neyman-Pearson  criterion-  It  is  described  in  some  detail  by  Burdic  (1984) 
and  many  other  authors;  in  this  report  only  a  summary  is  provided. 

When  a  source  is  not  present,  the  received  signal  power  ao  (t)  consists 
entirely  of  noise  n  (t)  : 

ao  (t)  =n(t)-  (28) 

When  the  source  is  present,  the  received  power  ao  (t)  contains  the  source 
signal  s  (t)  as  well  as  noise: 


ai  (t)  =  s(t)  +  n  (t)  - 


(29) 
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The  received  signal  power  in  each  of  these  cases  has  a  distinct  probability 
density  function  (pdf).  We  indicate  the  pdf's  for  the  two  cases  as  po  (a)  arid 
Pi  (a),  respectively 

Suppose  we  select  a  threshold  7,  below  which  it  is  assumed  that  no  source 
is  present,  and  above  which  it  is  assumed  that  a  source  is  present.  The 
probability  of  a  false  alarm  is  then 

roo 

Pfa  =  Po  (a)  da.  (30) 

J7 

The  probability  of  detection  is 

roo 

Pd=  Pi  (a)  da,  (31) 

“'7 

and  the  probability  of  false  dismissal  is 

Pfd  =  [  pi{a)  da  =  l-  Pd-  (32) 

Jo 

The  Neyman-Pearson  criterion  consists  of  fixing  Pfa  to  some  desired 
value,  and  then  determining  the  threshold  7  corresponding  to  this  value 
of  Pfa,  using  equation  (30).  This  procedure  obviously  requires  a  model  for 
Po  (o)-  Once  the  value  for  7  has  been  obtained,  we  can  compute  Pd  from  a 
model  for  pi  (a)  using  equation  (31). 

In  general  it  is  quite  a  difficult  problem  to  determine  the  pdf's  for  signals 
propagated  through  the  atmosphere.  Let  us  assume  initially  that  the 
signals  are  fully  saturated.  This  means  that  atmospheric  turbulence  has 
broken  up  the  deterministic  signal  into  a  large  number  of  statistically 
independent  contributions  having  random  phases.  The  statistics  for  such 
strongly  scattered  signals  are  well  known.  Specifically,  if  the  signal  is 
regarded  as  complex  (as  is  conventional  in  acoustics  and  many  other  fields 
involving  wave  propagation),  the  real  and  imaginary  parts  are  zero-mean, 
Gaussian  random  variables  (rv's)  having  the  same  variance  cj^/2 
(Tatarskii,  1971).  The  complex  signal  corresponding  to  s  (t),  for  example,  is 
written 

s  (t)  =  \Js{t)  exp  [i  {ujt  -h  (j)  (i))]  =  sr  (t)  -h  isj  (t)  ,  (33) 

where  s  is  the  signal  magnitude  squared,  f  is  the  phase,  u)  =  27r/  is  the 
angular  frequency,  sr  is  the  real  part,  and  sj  is  the  imaginary  part.  Note 
that  s  equals  the  sum  of  the  squares  of  two  independent  Gaussian  rv's: 

s(t)  =  s  (t)  s*  (t)  =  si  (t)  +  sj  (t)  . 
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Figure  7.  Probability 
density  function  for  a 
chi-squared  random 
variable.  Mean  is  set  to 
one,  and  parameter  is 
number  of  degrees  of 
freedom. 


(The  asterisk  indicates  the  complex  conjugate.)  It  can  be  shown  that  s  in 
this  case  has  an  exponential  pdf: 


where  H  (x)  is  the  Heaviside  function  (0  for  x  <  0  and  1  otherwise),  and 
cj^  is  the  variance  of  s.  It  can  be  shown  that  (s)  =  a^.  ABFA  uses  this 
relationship  to  determine 

For  the  pdf  of  the  noise,  ABFA  uses  a  chi-squared  distribution.  A 
chi-squared  pdf  of  order  u  corresponds  to  the  sum  of  u  independent, 
zero-mean  Gaussian  random  variables,  each  having  the  same  variance 
It  is  given  by 


(35) 


It  can  be  shown  that  the  mean  value  of  a  for  the  chi-squared  distribution  is 
(a)  =  ua^.  Hence  if  the  mean  signal  level  and  u  are  specified,  the  variance 
becomes  fixed.  The  chi-squared  pdf  is  plotted  for  several  values  of  u  in 
figure  7,  with  (o)  =  1.  Note  that  as  u  is  increased,  the  pdf  becomes  a  delta 
function  centered  on  the  mean.  The  limit  v  — >  oo  therefore  corresponds  to 
a  constant  noise  background. 

Note  that  when  v  =  2,  the  chi-squared  pdf  reduces  to  an  exponential  pdf. 
Physically,  one  might  think  of  v/2  roughly  as  the  number  of  independent. 


Figure  7.  Probability 
density  function  for  a 
chi-squared  random 
variable.  Mean  is  set  to 
one,  and  parameter  is 
number  of  degrees  of 
freedom. 
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complex  signals  contributing  to  the  total  noise.  If  there  are  very  few  noise 
sources  (small  u),  the  mean  intensity  can  undergo  dramatic  fluctuations.  If 
there  are  many  noise  sources  (large  u),  the  noise  background  is  nearly 
constant.  (Think  of  the  sound  created  by  a  few  crickets  chirping,  as 
opposed  to  a  field  full  of  crickets.)  It  should  be  stressed,  though,  that  one 
should  not  be  too  literal  in  interpreting  iv/2  as  the  number  of  independent 
noise  sources.  The  assumptions  inherent  to  the  chi-squared  pdf  (strong 
scattering,  independent  noise  sources  of  equal  variance)  are  unlikely  to  be 
satisfied  in  a  given  situation.  Therefore  it  is  more  reasonable  to  view 
equation  (35)  as  a  qualitative  model  whose  parameters  need  to  be 
determined  empirically. 

Let  us  now  determine  the  probability  of  false  alarm  using  the  chi-squared 
pdf.  The  integral  we  need  to  solve  is 


Pfa  = 


p  {a/aiyP  ^ 
i  a22^/2r(v/2)""P 


_  M 

2alJ 


da. 


By  defining  a'  =  al2a\,  we  have 

1 


P!a  = 


r  (i^/2)  J'y/2al 


[  (a'yP  ^  exp  (—o')  da. 
U/2ul 


The  integral  can  now  be  expressed  easily  in  terms  of  the  incomplete 
gamma  function,  defined  as  (Abramowitz  and  Stegun,  1965,  eq  (6.5.1)) 

P{a,x)  =  -^f  e~H^~^dt,  (36) 

r  (a)  To 


where  L  ()  is  the  ordinary  (complete)  gamma  function.  The  result  is 


Pfa  =  l-P 


2’2a^J- 


(37) 


There  is  no  way  to  solve  this  equation  for  7  in  closed  form.  However,  it  is 
fairly  easy  to  solve  using  numerical  methods  such  as  the  bisection 
technique. 

Next  let  us  determine  the  probability  of  detection.  Since  the  signal  in  this 
case  is  the  sum  of  two  rv's,  the  pdf  pi  (a)  is  found  by  convolution  of  the 
pdf's  of  the  two  rv's  (Bendat  and  Piersol,  1986): 

/OO 

po  (z)p(a  -  z)  dz. 

-OO 


Substituting  this  integral  into  equation  (31)  and  changing  the  order  of 
integration,  we  find 


/OO  r  poo 

Po  (z)  /  p{a-  z)  da 

-OO  {.Jo 


dz. 
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Using  equation  (34),  we  evaluate  the  inner  integral  as 


z)  da 


exp  [-  (7  -  z)  /cr^]  ,  z  Kj, 
1,  z>-f. 


When  this  result  is  substituted  back  into  the  outer  integral,  we  have 


Pd  = 


2alY{yl2) 


exp 


-i.\ 


1//2-1 


a'^J  Jo 


exp 


— 2; 


(tV2 


dz 


+ 


(a) 


i^/2-l 


exp 


-a) 


dz 


Setting  z”  =  zj 2a‘l  in  the  first  integral,  where 


2  ^ 

"  -  2cj2 


(38) 


and  z’  =  z/ 20-2  in  the  second  integral,  we  have 


Pd  = 


exp 


t2 


r(v/2) 

+  /  izY^‘^~^  eyi^)  i—z)  dz' 

Jznal 


^  /  (-2  )  exp(-z  )  dz 

\(ji  Jo 


''yl2o 

The  integrals  are  easily  evaluated  in  terms  of  the  incomplete  gamma 
function  (eq  (36)),  with  result 


Pd  =  l-P 


=  P^,  +  exp(-^l  ^ 


+  exp|-4U4V'^ 


2'2al) 


v/2 


P 


2'2al)- 


(39) 


Note  that  the  probability  of  detection  is  always  greater  than  the 
probability  of  false  alarm.  This  is  because  the  power  in  the  signal  is  always 
higher,  on  average,  when  the  source  is  present. 

Some  example  curves  for  Po,  as  a  function  of  P/a,  are  shown  in  figure  8. 1 
created  these  curves  by  determining  the  threshold  7  from  P/a,  using 
equation  (37).  The  mean  signal  level  (s)  was  set  to  1.  The  mean  noise  level 
(n)  then  follows  from  the  SNR  in  decibels  according  to 

SNR  =  101ogj4.  (40) 

(n) 
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Figure  8.  Probability  of 
detection  Po  as  a  function 
of  probability  for  false 
alarm  Pfa  for  a  single 
sample  of  received  signal. 
An  exponential  pdf  was 
used  for  source  level 
fluctuations,  whereas  a 
chi-squared  pdf  was  used 
for  backgrond  noise. 
Signal-to-noise  ratio  (SNR) 
and  degrees  of  freedom  in 
chi-squared  pdf  are 
parameters. 


Last,  the  variances  follow  from  the  relations  (t2  =  (s)  ,  and  =  (n)  tv- 

The  curves  in  figure  8  show  that  Pd  increases  with  increasing  Pfa,  as  one 
would  expect.  It  is  also  clear  that  detection  is  much  better  when  the  SNR  is 
high.  Increasing  v  favors  detection  somewhat,  particularly  at  low  SNR.  It 
should  be  pointed  out  that  the  highest  values  of  Pfa  on  the  plot,  0.1, 
would  be  much  too  high  in  most  practical  situations.  For  example,  if  an 
effort  is  made  every  1  min  to  detect  a  source,  false  alarms  would  occur 
every  10  min  on  the  average.  It  would  be  impractical  to  investigate 
possible  detections  so  frequently.  Hence  false  alarm  rates  of,  say,  10“^  or 
less  are  much  more  realistic. 

7.3  Probability  of  Detection:  Multiple  Receivers 

The  analysis  in  the  previous  section  assumed  that  just  one  receiver 
(sensor)  is  used  to  detect  the  source.  But  what  happens  if  data  from 
multiple  receivers  are  merged? 

It  is  actually  easiest  to  address  the  multiple-receiver  problem  in  terms  of 
probabilities  of  false  dismissals.  Suppose  there  are  N  receivers.  Let  us 
indicate  the  probability  of  false  dismissal  for  one  of  the  receivers  i  as 
Pfd,i  =  1  —  Pd^i-  Since  the  probability  of  occurrence  for  multiple  events  is 
found  by  multiplication  of  the  probabilities  for  the  individual  events,  the 
probability  of  false  dismissal  by  all  the  receivers  is 

N 

Pfd  =  X{Pfd,.  (41) 

i=l 
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(42) 


Hence  the  probability  of  detection  by  one  of  the  receivers  is 

N 

-pd = 1 — Pfd = 1  “  n  (1  “  P-D,*)  • 

i=l 

In  addition  to  the  possibility  of  multiple  receivers,  we  need  to  account  for 
the  fact  that  a  single  receiver  may  be  used  to  collect  multiple  samples.  This 
problem  is  quite  similar  to  the  multiple-receiver  problem.  Suppose  that  a 
total  of  Mi  statistically  independent  samples  are  collected  and  analyzed  by 
receiver  i.  If  the  source  of  interest  is  not  present,  the  probability  of  not 
having  a  false  alarm  in  a  single  sample  is  (1  —  P/a,*)/  arid  the  probability  of 
having  no  false  alarms  at  all  is  (1  —  P/a,i)^T  Hence  the  probability  that 
there  will  be  any  false  alarms  among  the  Mi  samples  is 

Pfa,^  =  1  -  (1  -  P/a,*)"^*  .  (43) 

Similarly,  if  Mi  samples  are  collected  while  the  source  is  present,  the 
probability  that  the  source  will  be  detected  among  any  of  the  samples  is 

PD,^  =  1  -  (1  -  Pd,*)"^*  .  (44) 

In  a  given  scenario,  it  would  normally  be  the  trade-off  between  Pfa.^  i  and 
Pd,*  that  is  of  interest,  as  opposed  to  the  trade-off  between  Pfa,i  and  Pd,*. 
(Note  that  the  two  are  the  same,  however,  when  Mi  =  1.)  Hence  it  is 
appropriate  to  use  Pd,*,  rather  than  Pd,*,  in  equation  (42).  We  can 
determine  the  relationship  between  Pfa,i  and  Pd,*  for  a  specified  Mi  by 
first  solving  (43)  for  Pfa,i' 

P/a,*  =  l-(l-P/a,*)'^'^P  (45) 


This  value  for  Pfa,i  can  then  be  used  in  the  single-receiver  analysis 
described  in  the  previous  section.  The  resulting  value  for  Pd,*  is  then 
simply  substituted  into  equation  (44).  Figure  9  shows  the  results  of  such  a 
calculation,  for  several  values  of  M*  (1,  2, 4,  and  10).  Clearly,  detection  is 
enhanced  when  multiple  independent  samples  can  be  collected. 


The  number  of  independent  samples  M*  is  related  to  the  time  duration  % 
of  the  data  records  processed  by  the  sensing  algorithm.  The  quantity  %  is 
called  the  sensor  integration  time  by  AFBA.  ABFA  estimates  M*  by  dividing 
Ti  by  the  turbulence  time  scale.  For  propagation  near  the  ground,  the 
appropriate  turbulence  time  scale  is  z/rt*.*  Hence 


M*  ~  max 


(46) 


*The  equation  for  the  turbulence  time  scale  is  a  very  rough  estimate,  and  is  valid  only 
when  wrnd-shear-generated  turbulence  dominates  the  scattering.  This  equation  results 
from  dimensional  analysis,  with  a  being  the  relevant  length  scale  in  a  shear  layer,  and  u* 
(the  friction  velocity)  being  the  relevant  velocity  scale. 
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Figure  9.  Probability  of 
detection  Po  as  a  function 
of  probability  for  a  false 
alarm  Pfa  for  multiple 
samples  of  received  signal. 
SNR  is  0  dB,  and  degrees 
of  freedom  in  background 
noise  is  u  =  16.  M  is 
number  of  statistically 
independent  samples. 


The  reason  for  placing  a  lower  bound  of  1  on  the  value  of  Mi  is  that  any 
sample  by  the  sensor  produces  at  least  one  statistically  independent  data 
point.  When  the  sensor  listens  for  intervals  longer  than  the  turbulence 
time  scale,  more  independent  samples  are  produced. 


7.4  Probability  of  Detection:  Unsaturated  Signals 

In  section  7.2  it  was  assumed  that  the  signal  to  be  detected  was  fully 
saturated.  This  assumption  allowed  us  to  derive  an  equation  for  the 
probability  of  detection  in  closed  form.  More  generally,  signals  may  have 
varying  degrees  of  saturation.  The  pdf  for  an  unsaturated  signal  has  a 
Rice-Nakagami  probability  density  function  (Platte  et  al,  1979): 


p{s)  =  ^exp 


s  +  m?  \ 


lo 


H{s), 


(47) 


where  Iq  ()  is  the  modified  Bessel  function.  In  the  limit  m  ^  0,  equation 
(47)  reduces  to  (34). 

From  equation  (32)  and  the  fact  that  the  pdf  of  the  sum  of  two  rv's  is  given 
by  the  convolution  of  the  pdf's,  we  find 


Pd 


1 

1 


Pi (a)  da 


z)  p  {z)  dz 


da. 
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Changing  the  order  of  integration,  and  using  the  property  that  both  po  arid 
p  are  zero  when  their  arguments  are  less  than  zero,  we  obtain 


Pd  =  1- 


z)  da 


p  (z)  dz. 


If  the  pdf  for  the  noise  intensity  is  still  assumed  to  be  chi-square,  as 
specified  by  equation  (35),  the  inner  integral  can  be  expressed  as  an 
incomplete  gamma  function,  and  we  have 


This  equation  actually  is  valid  for  any  p  (s).  When  equation  (47)  is 
specifically  used  for  p  (s),  the  integration  must  be  done  numerically  This 
makes  calculation  of  the  detection  probability  for  unsaturated  signals 
more  computationally  intensive  than  for  fully  saturated  signals. 

The  parameter  m  can  be  determined  with  equation  (8.3.9)  from  Platte  et  al 
(1979),  after  the  equation  is  rescaled  by  the  signal  intensity: 


m  = 


(49) 


In  this  equation,  4>  is  called  the  strength  parameter  of  the  turbulence.  The 
strength  parameter  is  given  by 

=  (50) 

where  /cq  is  the  acoustic  wavenumber,  R  is  the  propagation  distance,  L  is 
the  integral  length  scale,  and  is  the  index-of-refraction  variance. 

ABFA  uses  the  atmospheric  model  (sect.  4)  to  estimate  the  strength 
parameter.  Once  m  has  been  determined,  the  signal  variance  parameter 
follows  from  the  relationship  {s)  =  +  m?,  which  can  be  proven  with 

equation  (47). 


Target  Direction-Finding  Accuracy 

Acoustic  sensor  arrays  determine  the  direction  of  a  source  by  estimating 
the  angle  at  which  the  acoustic  wavefronts  arrive  at  the  array.  The  method 
used  by  ABFA  for  determining  the  accuracy  of  acoustical  bearing 
estimates  is  discussed  in  detail  elsewhere  (Wilson,  1997a).  Here  only  a 
summary  of  the  method  is  given. 


It  can  be  shown  that  the  variance  of  the  angle-of-arrival  estimates  ip  about 
the  actual  value  ip  is  always  greater  than  or  equal  to  the  inverse  of  a 
quantity  known  as  the  Fisher  information  matrix  J : 

>  J-'W-  (51) 

(In  the  case  considered  here,  involving  only  a  single  source,  J  is  a  scalar.) 
The  right-hand  side  of  equation  (51)  is  called  the  Cramer-Rao  lower  bound 
(CRLB).  Song  and  Ritcey  (1996)  show,  for  signals  having  a  joint-Gaussian 
probability  distribution,  that 

JW  =  A/tr(R-i^R-i^),  (52) 

where  M  is  the  number  of  statistically  independent  samples  (see  sect.  7.2), 
and  Raa  is  the  cross-correlation  matrix  for  the  signals  received  at  each  of 
the  sensors,  defined  by 


Raa  {lp)  =  (a(V’,  t)a{lp,  t))  ,  (53) 

in  which  a(V',  t)  is  a  column  vector  whose  elements  are  the  time-varying 
acoustic  pressure  signals  at  each  of  the  N  sensors  in  the  array.  As  before, 
we  can  decompose  Sipip,  t)  into  contributions  from  the  source  of  interest 
and  noise: 

a(V’,  t)  =  s{ip,  t)  -\-  n{t).  (54) 

Assuming  that  the  source  signal  and  noise  are  uncorrelated,  the  received 
signal  cross-correlation  matrix  becomes  the  sum  of  R^s  {'ip)  = 

{s{ip,t)s{ip,t))  andRnn  =  (n(i)n(t)): 

Raa  (V’)  =  R55  (V')  +  Rnn-  (55) 

Assuming  furthermore  that  the  noise  at  the  individual  sensors  is  mutually 
uncorrelated  and  equal  in  variance,  we  have 

Raa  (ip)  =  Rss  (V’)  +  0-nI>  (56) 

where  is  the  noise  variance.  When  the  signals  and  noise  are  normalized 
so  that  s{ip,  t)  has  unit  variance,  becomes  a  noise-to-signal  variance 
ratio.  It  is  related  to  the  SNR  (in  decibels)  by 

SNR  = -lOlogu^.  (57) 

Following  Wilson  (1997a),  the  source  signal  correlation  matrix  is  written  in 
the  form 

R,,  (iP)  =  S{iP)Q  T, 
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where  S  (V^)  is  called  the  array  steering  matrix,  and  T  accounts  for  the 
imperfect  coherence  of  the  signal  across  the  array  due  to  turbulence.  The 
symbol  ©  indicates  the  Hadamard  matrix  product  (simple  element-by- 
element  multiplication).  The  steering  matrix  depends  on  the  phase  delays 
between  the  sensors;  the  element  in  the  mth  row  and  nth  column  is  the 
phase  delay  between  sensors  m  and  n  in  the  array.  This  delay  is  given  by 

Smn  (ip)  =  exp  [-ikdmn  COS  (V’  “  Omn)]  ,  (58) 


where 

^mn  ~  ~  Xm)  “I"  {Vn  ~  Uva) 

is  the  sensor  separation,  and  (xn,  Vn)  is  the  location  of  sensor  n; 


amn  =  arctan  [{yn  -  Vm)  /  {xn  -  x^)] 


is  the  angle  between  the  sensors. 

The  elements  of  the  matrix  describing  the  turbulence  effect  are  equal  to  the 
mutual  coherence  function  (MCF),  evaluated  at  the  sensor  separation: 

Tmn  =  T  (dmn)  ■  (59) 


The  MCF  is  discussed  in  section  6.6. 


7.6  Target  Location-Finding  Accuracy 

When  the  data  from  two  or  more  direction-finding  arrays  are  merged,  it  is 
possible  to  estimate  the  actual  location  of  a  source.  ABFA  calculates  the 
uncertainty  in  such  location  estimates  using  a  very  general,  inverse 
problem  formulation  (Aki  and  Richards,  1980;  Wilson  and  Thomson, 

1994). 

Suppose  that  the  angles-of-bearing  estimates  from  N  sensor  arrays  are 
used  to  estimate  the  location  of  the  source.  The  locations  of  the  N  arrays, 
{xn,  Vn),  n  =  1,2, . . .  ,N,  are  assumed  to  be  known  exactly.  The  location  of 
the  source  (target),  {x,  y),  is  assumed  to  be  known  to  within  some 
uncertainty  a  in  both  the  x-  and  y-directions.  (In  ABFA,  the  initial 
uncertainties  in  the  source  position  are  set  to  the  overall  size  of  the 
domain.  Note  that  the  notation  in  this  section  differs  somewhat  from  that 
in  the  previous  one.)  The  full  geometry  of  the  problem  is  illustrated  in 
figure  10. 

The  angle-of-bearing  estimate  produced  by  each  of  the  arrays  is  used  to 
reduce  the  initial  uncertainty.  Let  ipn  be  the  actual  bearing  of  the  source 
relative  to  each  receiving  array: 

tanV’n  =  — — —■  (60) 
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Figure  10.  Geometry  of 
inverse  problem  for 
location-finding  accuracy. 


assumed 


Now  let  us  define  ■0^  as  the  estimated  (reference)  direction  of  the  source, 
and  (x,  y)  its  estimated  location.  The  perturbations  of  the  estimated 
position  from  the  actual  one  are  defined  as 

=  V'n  -  X  =  X  -  X,  v'  =  V  -  V- 

In  the  context  of  inverse  theory,  the  "knowns"  (the  are  usually  called 
the  data.  Arranging  the  data  as  a  column  vector,  we  define 


V’i 


V’jv 


The  unknown  perturbations  x'  and  y'  to  be  determined  are  called  the 
models: 


m  = 


The  so-called  forward  problem  consists  of  determining  the  data  from  the 
models.  Equation  (60)  gives  the  solution  to  the  forward  problem.  The 
inverse  problem  consists  of  solving  for  the  models  in  terms  of  the  data. 
Although  solution  of  the  inverse  problem  is  far  from  trivial,  it  can  be  made 
much  more  tractable  if  we  first  linearize  the  forward  problem,  by  keeping 
only  the  initial  terms  in  the  multivariate  Taylor  series  expansions  for  the 
data.  That  is. 


V’n 


dx 


x 


dipr 


x,y 


dy 


^,y 
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The  forward  problem  can  be  written  in  matrix  notation  as 


Vx  (x  -  Xn) 


d  =  Gm  +  n, 


where  the  vector  n  represents  the  noise  (uncertainty)  in  the  data,  and 


G 


dipi/dx  dipi/dy 


dil^N/dx  d'lpN/dy 


From  the  forward  problem  formulation,  equation  (60),  we  find 


(61) 


dll^n 

dx 


dlpn 


dy 


x,y 


cos 

X-  Xn 
sin  cos  ipn 

y-Vn 


(62) 

(63) 


Let  us  now  consider  the  somewhat  more  complicated  problem  where  the 
source  is  moving  at  a  velocity  {vx,Vy).  The  bearing  to  the  source  is  then 
given  by  equation  (60),  but  with  x  and  y  replaced  by  the  location  where 
the  sound  was  emitted.  This  location  is  {x  —  Vxt,  y  —  Vyt),  where  t  is  the 
time  of  emission.  Hence 


tan  -ipn 


y  -Vyt-  yn 
X  —  Vxt  —  Xn 


(64) 


We  can  determine  the  time  of  emission  by  equating  cot  to  the  propagation 
distance: 

Cot  =  \J{x-  Vxt  -  Xnf  +  {y-  Vyt  -  ynf. 

Solving  this  equation  for  t,  one  finds 


+  Vy{y-  yn)  +  y [Vx  {x  -  Xn)  +Vy{y-  yn)f  +  (eg  -vl- 


(x  -  Xn)^  +  {y-  ynf 


Co-vt^-  v^ 


The  matrix  G  for  the  forward  problem  generalizes  to 


(65) 


G 


dipi/dx  dipi/dy  dipi/dvx  dipi/dvy 
dipN/dx  dipw/dy  d^N/dvx  dipw/dvy 


(66) 


where 


5V’n 


dx 

9V’n 


(^sin  +  Vy/coj  cos  ■0^ 


dy 

9V’n 


dv^ 

d'lpn 


x,y,Vx,Vy 


x,y,Vx,Vy 


x,y,Vx,Vy 


X  —  Xr. 


sin'0„  (cos^„  +  Vx/co 


y-Vn 


sin  ipj.^ 


dv,, 


Co 

cosipn 

Co 


,  and 


\x,y,Vx,Vy 

In  ABFA,  the  estimated  velocity  {vx,Vy)  is  set  by  the  user. 


(67) 

(68) 

(69) 

(70) 


If  G  were  an  invertible  matrix  and  the  noise  were  negligible,  we  could 
now  solve  the  inverse  problem.  Unfortunately,  G  is  generally  not 
invertible;  in  fact,  it  is  not  necessarily  even  a  square  matrix.  Furthermore, 
uncertainties  in  the  data  are  often  substantial.  Fortunately,  methods  for 
solving  inverse  problems  such  as  this  have  been  studied  for  decades  in 
such  areas  as  seismology,  ocean  acoustic  tomography,  and  medical 
tomography.  One  of  the  simplest  and  most  satisfactory  methods  for 
approximating  the  inverse  is  called  the  stochastic  inverse,  which  minimizes 
the  expected  mean  square  difference  between  the  model  estimates  and 
their  actual  values.  (See,  for  example,  Aki  and  Richards  (1980),  or  Wilson 
and  Thomson  (1994).)  The  stochastic  inverse  is  given  by 


ihs  =  Gg  M,  (71) 


where 

(72) 

and 

R-mm 

=  (mih) , 

(73) 

Rnn 

II 

(74) 

Rmd 

=  ^md^  =  RmmG, 

(75) 

Rc/d 

=  /dd\  =  GR,,nmG  + 

(76) 

The  four  equations  above  define  the  model  covariance  matrix,  the  noise 
covariance  matrix,  the  model-data  covariance  matrix,  and  the  data  covariance 
matrix. 


With  regard  to  our  problem  involving  acoustic  source  localization,  we 
may  assume  that  the  initial  uncertainties  in  the  x-  and  y-coordinates  of  the 
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source  are  independent  and  both  equal  to  a  prescribed  value  The  initial 
uncertainty  in  the  velocities  is  Hence 


R 


mra 


0-2  0  0  0 

0  0-2  0  0 

0  0  0-2  0 

0  0  0  cj2 


(77) 


The  "noise,"  or  data  uncertainty,  in  this  case  is  given  by  the  CRLBs  of  the 
angle-of-arrival  estimates,  equation  (51): 


R 


nn 


J-MV-i)  0 

0  J-i  (V^s) 


0 

0  J-^(V'v) 


(78) 


Our  interest  in  an  application  such  as  ABFA  is  not  in  the  estimated  source 
positions  themselves,  but  rather  in  the  statistical  uncertainties  of  the 
estimates.  The  uncertainty  is  given  by  the  error  covariance  matrix 


Ree  =  (ee)  , 


(79) 


where  e  =  m  —  ih.  It  can  be  shown  (Wilson  and  Thomson,  1994)  that 

Ree  =  (W  -  I)  R,„,„(w  -  l)  +  G;iR„n^S  (80) 

where  W  =  is  called  the  resolution  matrix.  In  our  case,  Ree  is  a  4  x  4 

matrix.  The  upper  left  element  of  Ree  is  the  expected  variance  of  the  x 
estimates  about  their  actual  values,  whereas  the  lower  right  element 
corresponds  to  the  variance  in  y. 

To  summarize,  we  have  the  following  procedure  for  calculating  the 
uncertainties  in  the  source  location  estimates: 


1.  A  source  position  is  hypothesized.  The  uncertainty  of  the  source 
actually  residing  at  this  position  is  set  to  the  size  of  the  simulated 
domain.  Based  on  this  information,  the  model  covariance  matrix, 
equation  (  77),  is  determined. 

2.  Equations  (68)  and  (69)  are  used  to  determine  the  linearized  forward 
problem  solution  G,  equation  (61). 

3.  The  propagation  and  noise  models  are  run  to  calculate  the  SNR  at 
each  array  position,  for  a  source  at  the  hypothetical  position. 
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4.  The  SNR  is  used  to  calculate  the  CRLBs  for  the  direction-finding 
performance  of  each  array.  The  noise  covariance  matrix,  equation 
(78),  follows  from  the  CRLBs. 

5.  Now  that  Rmm/  Rnn/  and  G  are  known,  Hmd  and  follow 
immediately  from  equations  (75)  and  (76),  respectively.  The 
stochastic  inverse  is  then  determined  from  equation  (72). 

6.  The  reduced  uncertainty  in  the  source  position  (the  reduction  being 
due  to  the  information  provided  by  the  arrays)  is  given  by  equation 
(80). 
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8.  Example  Calculations 


8.1  Transmission  Loss 

As  discussed  in  section  6,  the  transmission  loss  (TL)  indicates  the  energy 
lost  by  the  sound  wave  as  it  propagates  from  the  source  to  the  receiver. 

The  TL  depends  strongly  upon  environmental  influences  such  as 
atmospheric  winds  and  temperature,  the  absorptive  properties  of  the 
ground,  and  the  topography  Some  example  calculations  are  provided  in 
this  section  to  demonstrate  these  environmental  influences. 

Figure  11  shows  a  TL  calculation  for  a  20-Hz  source.  This  frequency  is 
close  to  the  fundamental  frequency  produced  by  many  helicopters.  The 
height  of  the  source  is  10  m  above  ground  level.  The  height  of  the  sensor  is 
0.1  m,  and  its  horizontal  position  is  varied  across  the  simulation  domain. 
The  calculations  were  made  based  on  the  impedance  ground  model 
(sect.  6.2),  with  the  ground  type  being  long  grass.  Since  the  impedance 
ground  model  neglects  atmospheric  and  topographic  effects  on  the  sound 
levels,  the  sound  levels  decay  smoothly  and  omnidirectionally  away  from 
the  source. 

The  TL  calculation  shown  in  figure  12  is  essentially  a  repeat  of  figure  11, 
except  that  the  hybrid  wedge /impedance  ground  propagation  model 

Figure  11.  Transmission 
loss  calculation  for  a  20-Hz 
source  at  10-m  height 
(hovering  helicopter)  over 
grass,  using  plain 
impedance  ground  model. 

Location  of  helicopter  is 
marked  by  x . 

50 

¥ 

E  0 

D) 

'Q) 

JT 

-50 
0 


northing  (km) 


38 


(sect.  6.5)  was  used.  As  a  result,  topographic  effects  (but  not  atmospheric 
effects)  are  incorporated  into  this  calculation.  The  topography  was 
synthesized  with  the  von  Karman  model,  equation  (7),  for  L  =  2500  m, 
cj  =  20  m,  and  v  =  1.  Enhancement  of  the  sound  levels  is  evident  on  the 
sides  of  the  hills  facing  the  source,  whereas  acoustic  "shadows"  form  on 
hillsides  shielded  from  the  source. 

Figure  13,  like  figure  12,  uses  the  hybrid  wedge/impedance  ground 
model.  However,  the  surface  was  changed  from  long  grass  to  ice.  Hard 
materials  such  as  ice  and  rock  absorb  very  little  sound  energy.  As  a  result, 
the  TL  at  long  distances  from  the  source  is  much  reduced. 


Figure  12.  Same  as  figure 
11,  except  that  hybrid 
wedge /ground  impedance 
model  was  used.  Hence 
topographic  effects  on 
sound  field  are  included. 
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Figure  13.  Same  as  figure 
12,  except  that  groimd 
type  was  changed  from 
grass  to  ice. 
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The  next  set  of  TL  calculations,  figures  14  to  16,  are  again  for  a  20-Hz 
source  (helicopter)  at  10  m  over  grassy  ground,  except  that  atmospheric 
refraction  effects  are  incorporated  into  the  calculations  via  the  hybrid 
wedge/FFP  model  (sect.  6.5).  The  figures  show  calculations  for  three 
idealized  meteorological  conditions:  mostly  sunny  and  calm  (table  1), 
overcast  and  windy  (table  1),  and  a  deep  inversion  characteristic  of  the 
late  night  or  early  morning  under  clear  skies  (Noble,  1994).  For  the  mostly 
surmy  and  calm  conditions  (fig.  14),  there  is  a  sharp  increase  in 
temperature  near  the  surface.  As  a  result,  sound  energy  is  refracted 
upward,  resulting  in  a  large  TL  at  distances  greater  than  several 
kilometers. 


Figure  14.  Same  as  figure 
12,  except  that  hybrid 
wedge/FFP  propagation 
model  was  used.  Hence 
both  atmospheric  and 
topographic  effects  on 
sound  field  are  included. 
Atmospheric  conditions 
are  characteristic  of  a 
sunny,  calm  day. 
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Figure  15.  Same  as  figure 
14,  except  that  overcast, 
windy  atmospheric 
conditions  are  used.  Wind 
is  blowing  to  north  (from 
left  to  right,  and  slightly 
down,  in  figure). 
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Figure  16.  Same  as  figure 
14,  except  that  deep 
inversion  atmospheric 
conditions  are  used. 
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The  overcast,  windy  case  (fig.  15)  features  a  strong  wind  blowing  from  the 
north.  The  wind  speed  is  about  6.2  m/ s  at  a  height  of  10  m.  Because  of  this 
strong  wind,  the  TL  is  much  less  downwind  from  the  source  than  upwind. 
The  deep  inversion  case  (fig.  16)  has  a  temperature  gradient  of  +15  °C/km 
below  300  m,  and  —8  °C/km  above.  No  wind  is  included  in  this  case.  The 
positive  temperature  gradient  below  300  m  results  in  pronounced 
trapping  of  sound  energy  near  the  ground.  As  a  result,  low  TL  is  observed 
far  from  the  source.  "Skip  zones,"  or  isolated  ranges  of  diminished  sound 
levels,  are  also  evident.  The  skip  zones  result  from  combined 
refraction/ ground-reflection  effects. 


8.2  Signal-to-Noise  Ratio 

SNR  calculations  (sect.  7.1)  incorporate  the  source  level,  TL,  and 
environmental  noise.  Figures  17  to  19  show  SNR  calculations  for  the  same 
three  atmospheric  conditions  used  in  figures  14  to  16.  The  source  level  at 
1  m  in  each  case  is  170  dB,  which  is  characteristic  of  a  loud  source  such  as 
a  helicopter.  The  bandwidth  of  the  source  is  1  Hz.  The  background  noise  is 
set  to  30  dB  per  1/3  octave  band,  which  is  a  level  characteristic  of  a  quiet, 
natural  setting.  As  a  result,  most  of  the  noise  is  caused  by  turbulent 
pressure  fluctuations  on  the  microphone,  or  "wind  noise." 

Both  the  sunny,  calm  case  and  the  overcast,  windy  case  (fig.  17, 18)  have 
SNRs  favorable  to  detection  (above  0  dB)  within  a  radius  5  to  10  km  from 
the  source.  The  region  of  favorable  SNR  for  these  cases  is  small  not  only 
because  of  the  high  TL,  but  also  because  both  cases  feature  strong 
turbulence  and  hence  high  levels  of  wind  noise. 
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Figure  17.  SNR  calculation 
for  a  20-Hz, 

1-Hz -bandwidth  source  at 
10-m  height  (hovering 
helicopter)  over  grass, 
using  hybrid  wedge/FFP 
model.  Atmospheric 
profiles  are  characteristic 
of  sunny,  cahn  conditions. 
Noise  background  is 
characteristic  of  a  quiet, 
natural  setting.  Location  of 
helicopter  is  marked  by  x . 


Figure  18.  Same  as  figure 
17,  except  that  overcast, 
windy  atmospheric 
conditions  are  used. 
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The  deep  inversion  case  (fig.  19)  has  no  turbulence,  and  as  a  result  the 
SNR  is  favorable  for  distances  very  far  from  the  source.  (In  actuality, 
temperature  inversion  conditions  often  have  low-intensity  or  intermittent 
episodes  of  turbulence.  Because  turbulence  in  such  statically  stable 
conditions  is  poorly  understood  at  present,  ABFA  makes  the  simplification 
of  setting  the  turbulence  intensity  to  zero.  The  user  can  change  the 
turbulence  intensity  to  a  nonzero  value  if  desired.) 
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Figure  19.  Same  as  figure 
17,  except  that  deep 
inversion  atmospheric 
conditions  are  used. 
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8.3  Probability  of  Detection 

Probabilities  of  detection  depend  strongly  upon  the  SNR.  If  the  SNR  is 
well  above  0  dB,  the  source  can  generally  be  detected.  But,  as  discussed  in 
section  7.2,  other  factors  related  to  the  sensor  and  propagation 
environment  also  play  a  role  in  detection.  Among  these  are  the  integration 
time  of  the  sensor  (how  long  it  listens  for  the  source  signal),  and  the 
characteristic  time  of  the  signal  fluctuations. 

Figure  20  shows  the  probability  of  detection  for  the  same  helicopter-like 
source  considered  previously.  Sunny,  calm  atmospheric  conditions  were 
used  in  the  calculation.  Unlike  the  previous  figures,  where  the  source 
position  was  fixed  and  the  sensor  position  varied,  in  figure  20  the  sensor 
position  is  fixed  and  the  source  position  varies.  The  sensor  is  positioned  at 
a  height  of  0.1  m  above  ground  level,  on  an  elevated  region  near  the  center 
of  the  simulation  domain.  Because  of  the  elevated  position  of  the  sensor, 
terrain  effects  on  detection  are  rather  small.  The  sensor  integration  time 
used  in  the  calculation  was  T  =  10  s.  The  height  of  the  propagation  path 
was  taken  as  5.05,  which  is  the  average  of  the  source  and  sensor  heights. 
Since  the  friction  velocity  for  the  sunny,  calm  case  is  u*  =  0.1  m/ s,  the 
turbulence  time  scale  is  therefore  0.1/5.05  =  50.5,  and  the  number  of 
independent  samples  determined  from  equation  (46)  is 
M  =  max  (1,10/50.5)  =  1. 

Figure  21  is  similar  to  20,  except  that  the  sensor  integration  time  has  been 
increased  to  T  =  10  min.  This  larger  value  for  T  causes  M  to  increase  to  12, 
and  the  area  of  high  probability  of  detection  increases  as  a  result.  Figure  22 


43 


Figure  20.  Probability  of 
detecting  a  helicopter 
hovering  over  grass,  for 
sunny,  calm  atmospheric 
conditions.  Noise 
backgroimd  is 
characteristic  of  a  quief, 
natural  setting.  Sensor 
integration  time  is  10  s. 
Sensor  location  (fixed)  is 
marked  by  o. 
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Figure  21.  Same  as  figure 
20,  except  that  sensor 
integration  time  was 
increased  10  min. 
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is  also  similar  to  20,  except  that  sensors  at  two  locations,  on  opposite  sides 
of  a  hill,  have  been  included  in  the  calculation.  Together  the  two  sensors 
provide  good  detection  probability  over  most  of  the  simulation  domain, 
except  in  the  northwest  corner  and  a  valley  along  the  eastern  edge.  This 
figure  demonstrates  the  utility  of  ABFA  for  determining  good  sensor 
locations. 
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Figure  22.  Same  as  figure 
20,  except  that  two 
independent  sensors  are 
deployed. 


Figures  23  and  24  are  similar  to  figure  20,  except  that  different  weather 
conditions  were  used.  Specifically,  overcast  and  windy  conditions  were 
used  for  figure  23,  and  the  deep  inversion  profiles  for  figure  24.  The  sensor 
integration  time  for  both  figures  was  10  s,  which  resulted  in  M  =  1  in  each 
instance.  The  region  of  high  probability  of  detection  has  a  dramatic 
asymmetry  in  figure  23  due  to  the  strong  wind  effect.  Conditions  for 
detection  are  much  improved  when  the  source  is  upwind  from  the  sensor. 
The  probability  of  detection  is  high  everywhere  in  the  simulated  domain 
for  the  deep  inversion  case  (fig.  24),  as  a  result  of  the  high  SNR. 

8.4  Target  Location-Finding  Accuracy 

As  discussed  in  sections  7.5  and  7.6,  an  array  of  acoustic  sensors  can  be 
used  to  determine  the  bearing  of  a  target  (source).  Data  from  multiple 
arrays  can  then  be  combined  to  determine  the  target  location.  Some 
examples  of  the  accuracy  of  target  locations  that  can  be  obtained  from 
multiple  arrays  are  shown  in  figures  25  to  28.  In  each  figure,  six  separate 
sensor  arrays  have  been  deployed  within  the  simulation  domain.  Each 
array  consists  of  six  microphones,  equally  spaced  along  a  circle  of  radius 
1  m.  Calculations  are  presented  for  the  same  three  weather  conditions 
used  previously:  sunny  and  calm  (fig.  25  and  26),  windy  and  overcast 
(fig.  27),  and  deep  inversion  (fig.  28).  The  target  is  the  same  as  before:  a 
20-Hz  source  similar  to  a  hovering  helicopter. 


45 


1 


Figure  23.  Same  as  figure 
20,  except  that  overcast, 
windy  atmospheric 
conditions  are  used. 
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Figure  24.  Same  as  figure 
20,  except  that  deep 
inversion  atmospheric 
conditions  are  used. 
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The  target-location-finding  accuracy  is  clearly  best  for  the  deep  inversion 
conditions.  According  to  the  calculations,  the  location  of  the  target  can  be 
determined  to  an  accuracy  better  than  a  few  hundred  meters  over  much  of 
the  simulation  domain.  Although  the  location-finding  accuracy  for  the 
cases  with  high  turbulence  (sunny  and  calm,  and  overcast  and  windy)  is 
not  as  satisfactory,  useful  location  information  can  still  be  obtained.  The 
target  location  can  be  determined  to  better  than  1  km  in  regions  where  the 
sensors  are  placed  less  than  about  5  km  apart.  The  main  advantages  of 
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Figure  25.  Accuracy  for 
resolving  position  of  a 
target  (hovering 
helicopter),  under  sunny, 
calm  atmospheric 
conditions.  Sensors, 
marked  by  o's,  each  consist 
of  an  array  of  six 
individual  microphones. 
Sensors  are  positioned  at 
high  elevation  points. 
Noise  background  is 
characteristic  of  a  quiet, 
natural  setting. 
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Figure  26.  Same  as  figure 
25,  except  that  sensors  are 
positioned  along  a 
semicircular  valley. 


using  acoustic  sensors  instead  of  radar  in  a  scenario  such  as  this  would  be 
low  cost,  insensitivity  to  ground  clutter,  and  passive  operation  (i.e.,  the 
sensors  themselves  do  not  emit  a  detectable  signal). 

The  two  figures  corresponding  to  sunny,  calm  conditions  (25  and  26) 
contrast  location-finding  accuracy  from  two  different  sensor  layouts.  In 
figure  25,  the  arrays  are  closer  to  each  other,  and  positioned  on  terrain 
features  having  high  elevation.  In  figure  26,  the  arrays  are  placed  in  a 
low-elevation  region  having  the  appearence  of  a  long,  semicircular  valley 


47 


Figure  27.  Same  as  figure 
25,  except  that  overcast, 
windy  atmospheric 
conditions  are  used. 
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Figure  28.  Same  as  figure 
25,  except  that  deep 
inversion  atmospheric 
conditions  are  used. 


Although  the  spatial  coverage  of  the  different  sensor  configurations  has  a 
complicated  appearance,  it  can  be  understood  from  a  few  basic 
observations.  First,  a  single,  isolated  array  is  of  no  help  in  determining  the 
location  of  a  source;  it  can  only  determine  bearing.  Second,  when  an 
isolated  pair  of  arrays  occurs,  the  location  of  a  target  cannot  be 
determined  along  the  axis  of  the  pair,  since  no  triangulation  is  obtained 
from  the  directional  information  produced  by  each  array.  Therefore,  only 
when  three  or  more  arrays  are  close  enough  to  a  target  to  yield  accurate 
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directional  information  can  the  target  location  be  accurately  determined. 
For  the  layout  in  figure  26,  there  is  no  combination  of  three  or  more  sensor 
arrays  packed  closely  enough  to  yield  a  connected  region  of  good 
coverage.  On  the  other  hand,  the  layout  in  figure  25  features  five  arrays 
that  are  close  enough  to  provide  a  good,  connected  region  of  coverage.  A 
sixth  array  (positioned  at  a  northing  of  about  5  km,  and  easting  3  km)  is 
too  far  removed  from  the  other  five  to  be  useful  in  the  location  finding. 


49 


9.  Concluding  Remarks  and  Topics  for  Future  Research 


ABFA  synthesizes  many  recent  advances  in  computational  acoustics, 
software  capabilities,  and  information  theory  as  applied  to  acoustic  arrays. 
The  result  is  a  very  flexible  tool  for  evaluating  the  performance  of  acoustic 
sensor  configurations  in  various  environments.  Although  ABFA  is  a 
prototype,  and  perhaps  not  yet  sufficiently  refined  for  the  general  user,  it 
incorporates  many  ideas  that  can  serve  as  useful  starting  points  for  future 
development. 

The  prediction  methods  incorporated  into  ABFA  are  generally  state  of  the 
art.  Still,  there  remains  much  need  for  critical  examination  and  evaluation 
of  some  of  the  underlying  equations  and  techniques  used  for  propagation 
modeling.  In  all  likelihood,  there  are  scenarios  where  ABFA's  prediction 
methods  break  down.  Because  research  in  atmospheric  acoustics  has 
received  rather  modest  funding  during  the  past  several  decades,  and  the 
driving  application  has  often  been  noise  control  rather  than  detection, 
there  remain  many  scientific  issues  related  to  detection  and  signal 
coherence  that  need  to  be  addressed  by  future  research.  Some  of  these 
follow: 

1.  Calculation  of  mean  sound  levels  in  the  presence  of  turbulent  scattering.  It 
has  been  well  documented  that  turbulent  scattering  substantially 
raises  sound  levels  in  refractive  and  topographic  shadow  regions 
(Gilbert,  Raspet,  and  Di,  1990).  At  present,  methods  for 
incorporating  realistic  turbulent  scattering  into  acoustical 
calculations  are  quite  slow  and  impractical  for  use  in  a  tool  such  as 
ABFA.  Hence  efficiency  improvements  in  scattering  calculations  are 
much  needed. 

2.  Improvements  to  the  scattered  signal  statistics.  The  current  version  of 
ABFA  does  not  properly  account  for  the  statistics  of  the  scattered 
signal  in  the  so-called  "geometrical  acoustics"  regime,  which  is 
characterized  by  strong  phase  fluctuations  and  weak  amplitude 
fluctuations.  This  situation  occurs  most  often  for  low  frequencies 
and  short-to-moderate  propagation  distances.  Turbulent  intermittency 
(the  tendency  of  turbulence  to  occur  in  bursts  of  activity)  can  also 
significantly  affect  detection  (Wilson,  1996)  and  is  unaccounted  for  at 
present. 
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3.  Improved  calculation  of  the  signal  coherence.  The  signal  coherence 
model  in  ABFA  assumes  straight-line  propagation.  As  a  result,  such 
potentially  important  effects  as  refraction  and  ground  reflections  are 
neglected  in  the  coherence  calculation.  Future  research  on  acoustic 
signal  coherence,  both  experimental  and  theoretical,  is  needed. 

4.  Decorrelation  time  of  the  acoustic  signals.  ABFA  estimates  the  time  scale 
for  decorrelation  of  the  acoustic  signal  in  a  rather  crude  fashion  from 
the  friction  velocity  and  height  of  the  propagation  path.  This  is 
necessary  because  so  little  information  is  available  at  present  on  the 
temporal  decorrelation  of  acoustic  signals. 

5.  Reliance  on  structure-function  parameters.  Structure-function 
parameters  are  convenient  and  accurate  for  describing  the  inertial 
subrange  (small-scale)  structure  of  turbulence.  However, 
propagating  acoustic  waves  are  known  also  to  be  significantly 
affected  by  energy  subrange  (large-scale)  turbulence  structure. 
Research  is  under  way  at  ART  to  address  this  issue  (Wilson,  1998a). 

6.  Improved  noise  statistics.  Very  little  is  known  about  the  statistical 
distributions  of  noise  in  various  environments,  although  this 
information  is  extremely  important  in  probability-of-detection 
calculations.  The  limitations  of  the  chi-squared  distribution  used  by 
ABFA  need  to  be  ascertained. 

7.  Wind  noise  models.  Wind  noise  plays  a  critical  role  in  determining 
detection  at  low  acoustic  frequencies.  The  wind  noise  model 
currently  used  by  ABFA  is  empirically  based,  and  the  limits  of  its 
applicability  are  not  well  known.  Somewhat  better,  but  still 
imperfect,  models  are  available.  It  would  be  valuable  to  incorporate 
one  of  these  other  models  into  ABFA,  or  to  develop  a  new,  superior 
formulation  for  the  wind  noise. 
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Abbreviations  and  Acronyms 


ABFA 

Acoustic  Battlefield  Aid 

ARL 

Army  Research  Laboratory 

CRLB 

Cramer-Rao  lower  bound 

FFP 

fast  field  program 

GUI 

graphical  user  interface 

MCF 

mutual  coherence  function 

MO 

Monin-Obukhov 

pdf 

probability  density  function 

rv 

random  variable 

SNR 

signal-to-noise  ratio 

TL 

transmission  loss 
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